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FAR-FIELD FEATURES OF THE KELVIN WAKE 


A. Barnell and F. Noblesse 


David W. Taylor Naval Ship Research and Development Center 
Bethesda, Maryland 20084 


ABSTRACT 


The classical Kelvin wake trailing a ship advancing 
at constant speed in calm water is studied. In particular, 
asymptotic expressions for determining the far-field wake, 
at large distances behind the ship, are investigated 
numerically. This analysis, in which surface tension and 
nonlinearities are ignored, indicates that the divergent 
waves of a surface ship generally are infinitely steep at the 
track of the ship, even though their amplitude vanishes 
there; this result is theoretically possible because the 
wavelength of the divergent waves vanishes at the track of 
the ship. Inasmuch as infinitely-steep water waves cannot 
exist in reality, the linear no-surface-tension analysis 
presented in this study suggests that no divergent waves 
can exist within a certain region in the vicinity of the 
track of the ship, and that the Kelvin wake contains three 
distinct regions: (i) an inner region where only transverse 
waves can exist, (li) an outer region where both transverse 
and divergent waves are present, and (ili) an intermediate 
region at the boundary between the inner and outer 
regions where short steep divergent waves, as well as 
transverse waves, can be found. Numerical results for a 
simple bow form show that the inner region is quite 
narrow, and that the wavelength of the divergent waves at 
the boundary of the inner region is of the order of | to 
20 cm, depending on the speed of the ship. These results 
appear to be consistent with the narrow V-wakes observed 
in some SAR images of ship wakes. The analysis 
presented in this study may thus provide a partial 
explanation for these observations. However, the small 
wavelengths mentioned above indicate that the divergent 
waves in the vicinity of the track of the ship are likely to 
be significantly altered by surface tension, which should 
therefore be included in a more realistic analysis. 


INTRODUCTION 
It has been observed, e.g. [1,2,3], that images of 
ship wakes taken by the SEASAT SAR (Synthetic 
Aperture Radar) and other airborne SAR systems 
sometimes reveal long, narrow wakes extending back 
several miles behind ships. According to Swanson [3], 
“there appears a narrow, dark band along the centerline 


of the ship wake, widening slightly with distance away 


from the ship. Often accompanying these dark regions, at 
a slightly larger angle, are bright lines, which show up 
especially well on X-band SAR images’’. The apparent 
included angle of this narrow wake is much smaller that 
that corresponding to the Kelvin cusp line. 

Several theoretical explanations of the features 
observed in SAR images of ship wakes have been 
proposed [3,4,5,6]. The proposed explanations include 
interactions between the cross currents created by a ship in 
its wake and surface gravit\ waves [3,6] and the occurence 
of a sharp peak in the amplitude of the divergent waves 
in the Kelvin wake for a ship form having a large flare 
angle [5]. However, there currently appears to be no 
convincing electrodynamic or hydrodynamic explanation of 
the observations. It thus seems necessary to investigate in 
detail each one of the proposed possible explanations in 
order to determine whether one, or possibly a combination 
of several, of them does in fact provide a satisfactory 
theoretical interpretation of the SAR images of ship 
wakes. 

The classical Kelvin wake trailing a ship advancing 
at constant speed in calm water is examined in this study. 
Viscosity, surface tension and nonlinearities are ignored in 
this preliminary analysis, and irrotational flow is assumed. 
The analysis and related numerical results for a simple 
ship bow form have led to several conclusions and 
recommendations for further studies which are presented 
at the end of the paper. Only two conclusions that may 
be related to some of the features observed in SAR 
images of ship wakes are discussed here. 

It was previously found by Scragg [5] that, for a 
ship bow form with a large flare angle, the zeroth-order 


slender-ship approximation to the far-field wave-amplitude 


function given in Noblesse [7] predicts a sharp peak in the 
value of the amplitude of the divergent waves at an angle 
from the track of the ship equal to approximately half the 
bow entrance angle. This finding of Scragg has been 
verified in this study, as may be seen from figures 15, 16 
and 18. Furthermore, the magnitude of the wave- 
amplitude peak has been found to increase very rapidly as 
the Froude number decreases below a certain threshold 
value. This peak thus appears to be a large-flare low- 
Froude-number feature. For the simple ship bow form 
considered in this study, the peak in the amplitude of the 
divergent waves in the Kelvin wake occurs along a line at 
6° from the track of the ship. It may be found from 
equation (30a) that the wavelength of the divergent waves 
along this line varies between 0.7 m and 4.5 m for a ship 
with speed varying between 10 knots and 25 knots, 
respectively. 

The other conclusion of this study that may be 
related to the narow V-wakes observed in some SAR 
images of ship wakes is the result that the divergent waves 
of a surface ship are infinitely steep at the track of the 
ship, even though their amplitude vanishes there. This 
result is theoretically possible because the wavelength of 
the divergent waves vanishes at the track of the ship. A 
similar result was previously obtained by Sharma [8] who 
found that the Michell thin-ship approximation for a thin 
and deep strut-like ship form predicted infinite slopes for 
the divergent waves at the track of the ship. Inasmuch as 
infinitely-steep water waves cannot exist in reality, the 
foregoing result suggests that no divergent waves can exist 
within a certain region in the vicinity of the track of the 
ship, and that the Kelvin wake contains three distinct 
regions: (i) an inner region where only transverse waves 
can exist, (ii) an outer region where both transverse and 
divergent waves are present, and (iii) an intermediate 
region at the boundary between the inner and outer 
regions where short steep divergent waves, as well as 
transverse waves, can be found. 

Numerical results for a simple bow form show that 
the ‘‘no-divergent-wave’’ inner region is quite narrow, as 
may be seen from figure 20 showing the Kelvin cusp line 
(angle ~ 19°28’), the line along which the amplitude of 
the divergent waves exhibits a peak (angle ~ 6°), and the 
three lines along which the steepness of the divergent 
waves is equal to 1/20, 1/15 and 1/7 (chain line close to 
the track of the ship). The latter three lines, along which 


the divergent waves are steep, lie much closer to the track 


of the ship than the line along which the steepness of the 
divergent waves exhibits a peak. Figures 21 and 22 show 
considerable variations among the several ‘‘steep-divergent- 
waves’’ lines that are represented in these figures. 
Nevertheless, these lines may be seen to correspond to 
values of Y/(—X) equal to about 10~2 to 2.1072. For 
ship speeds varying between 10 knots and 25 knots, the 
wavelength of the divergent waves corresponding to the 
foregoing values of Y/(— X) may be shown to vary 
between 0.7 cm and 4.2 cm for Y/(—X) = 10-2, and 
2.7 cm and 17 cm for Y/(—X) = 2.10~2. These 
wavelengths are consistent with the wavelengths of the 
radar pulses used in SAR imaging, so that the foregoing 
results may provide a partial hydrodynamic explanation 
for the narrow V-wakes observed in these images. 

The tentative nature of this explanation must 
however be stressed. Indeed, the foregoing results are 
based on an analysis in which surface tension and 
nonlinearities have been neglected. Inasmuch as this linear 
no-surface-tension analysis predicts extremely short and 
steep waves in the vicinity of the track of the ship, it is 
evident that both surface tension and nonlinear effects are 
liable to be significant. In particular, the short 
wavelengths found along the steep-divergent-wave lines 
determined in this study, and the brief description of the 
effects of surface tension upon the Kelvin wake given in 
Sharma [8], Lamb [9, pp. 468-470] and Wehausen and 
Laitone [10, pp. 636-637] indicate that the system of 
divergent waves in the vicinity of the track of the ship is 
likely to be profoundly affected by surface tension. 
Effects of surface tension upon the Kelvin wake will be 


investigated in a sequel to the present study. 


APPROACH 

This study considers the steady potential flow due 
to a ship advancing with constant speed in calm water of 
infinite depth and lateral extent. The far-field Kelvin 
wake, which is of primary interest here, may be 
conveniently analyzed in terms of the nondimensional far- 
field coordinates X = Xg/U2, velocity potential ¢ = og/U3 
and velocity vector V,¢ = Vy®/U, where g is the 
gravitational acceleration and U is the speed of advance of 
the ship, X and © represent the dimensional coordinates 
and velocity potential, respectively, and V, and Vy are the 
nondimensional and dimensional differential operators 
V, = (8/0x,0/dy,0/dz) and Vy = (0/0X,0/4Y,0/4Z). 


The mean free surface is taken as the plane z = 0, with 


the z axis pointing upwards, and the x axis is chosen in 
the ship centerplane and pointing towards the bow. The 
origin of the system of coordinates is placed within the 
ship. The Froude number is denoted by F = U/(gL)!2, 
where L is the length of the ship. 

Equation (32) in [7] yields the following expression 
for the velocity potential associated with the Kelvin wake 
behind the ship 


n¢(x) = mf, [E (GX) +E_(txIK(@® dt, (1) 


where E ,, (13x) is the exponential function 
E (5x) = exp[2(1 +t?) +i tyr + t?)!7), (2) 
and K(t) is the far-field wave-amplitude function, which 
depends on the hull shape and the Froude number. 
Assuming that differentiation under the integral sign is 
permitted in equation (1), we may obtain 

$, fe (E,+E_ (1 +t2)!/2 (Ba) 
H. IN| aie if, Sch 11 +12)!/2 (3b) 


The nondimensional elevation e = Eg/U? of the 


free surface at a sufficiently-large distance behind the ship, 
such that nonlinearities may be neglected, is given by 
e(x,y) = 49(x,y,0)/dx. (4) 
The slopes of the free surface in the directions parallel 
and perpendicular to the ship course then are 

Be(x,y)/Ax = 074(x,y,0)/Ox2 (4a) 
Be(x,y)/dy = 324(x,y,0)/Oxdy. (4b) 


If differentiation under the integral sign in equation (3a) is 


permitted, we have 


o(E +E 1+? (5a) 
Tt Pas = -imf | a io 2 
dxy 0 {E,-E_ (+t?) (5b) 
The vertical velocity $, is given by 
ee (6) 


In this study, we are mostly interested in the value 
of the several flow variables defined above at the mean 
free surface and behind the ship, so that we have z = 0 
and x < 0. Expression (2) for the exponential function E, 


then becomes 


E, = expfix(1# at)(1+t?)'”], (7) 
where a is defined as 
a = —y/x. (8) 


For a ship with port- and starboard-symmetry, as is 
considered here, the Kelvin wake is symmetric about the 
ship track y = 0. We may then restrict the analysis of the 
Kelvin wake to the domain y 2 0 and x < 0, and assume 
aZ20. 

Two difficulties associated with the foregoing 


approach for numerically determining the potential, its 


gradient, and the free-surface slopes in the Kelvin wake 
are readily apparent and should be noted here. A first 
numerical difficulty stems from the oscillations of the 
exponential function E., given by equation (7), which are 
very rapid for large values of |x|. We have x = Xg/U2 
= (X/L)/F2. For a typical value of the Froude number 
equal to 0.2, say, we thus have x = 25 X/L; we then 
have x = —250 at 10 ship lengths behind the ship, and 
much larger values of |x| must be considered at greater 
distances behind the ship or/and for smaller values of the 
Froude number. Even for the comparatively-moderate 
value of x equal to —50, figure 1 shows that the 
functions E, and E _ oscillate quite rapidly. More 
precisely, figure 1 depicts the real parts of the functions 
E , (t;x,@) and E_(t;x,a) for x = —50; a = .1, .2, 
1/23/2, .4 and .5; and for 0 <t <7 and 0 <t <3 on 
the left and right sides, respectively. Figure 1 also 
indicates that the behavior of the function E , (t;x,a) 


strongly depends on the value of a. 
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Fig. 1 — Real Parts of the Functions E , (t;x,a) and 
E_(tx,a) for x = —50 anda = .1, .2, 1/23/, .4 and .5 


A second, more basic, difficulty is associated with 
the differentiation under the integral sign which was used 
for obtaining expressions (3a,b) and (Sa,b) from 
expression (1). For a fully-submerged body, the far-field 
wave-amplitude function K(t) is exponentially-small as t > 
©, so that the operation of differentiating under the 
integral sign in expression (1) can be continued indefinitely 
in principle. Differentiation under the integral sign likewise 
is justified if z < 0. However, the operation must be 
justified in the limiting case z = 0 for a surface ship. 
Clearly, the operation may not be justified in principle, or 


feasible in practice, if the far-field wave-amplitude 


function K(t) does not vanish sufficiently rapidly as t > 
co. Precise information about the asymptotic behavior of 
the function K(t) as t > © is required in this respect. 
Figure 2 depicts the real and imaginary parts of the 
functions (1 +t2)K(t) and t(1 +t2)K(t) appearing in the 
integrands of the integrals (Sa) and (5b), respectively, for 
0< t< 16 and for a simple ship-bow shape which is 
considered further on in this study. Differentiation under 
the integral sign in expression (3a) is clearly not justified 


in the case corresponding to figure 2. 
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Fig. 2 — Real and Imaginary Parts of the Functions 
(1+t2)K(t) and t(1+t?)K(t) for a Simple Ship Bow Shape 
in the Zeroth-Order Slender-Ship Approximation 


ASYMPTOTIC EVALUATION OF THE 
KELVIN WAKE 


For large values of |x| and z = 0, analytical 
approximations to the integrals (1), (3a,b) and (Sa,b) can 
be obtained by taking advantage of the rapid oscillations 
of the exponential functions E , defined by equation (7). 


These functions may be expressed in the form 


E = explix0(t:a)], (9) 
where @(t;a) is defined as 
O(tsa) = (1—at)(1+t2)!/2, (10) 


The functions E , (t;x,a) and E_(t;x,a), where a > 0, 
correspond to the function E(t;x,a) with a > 0 and a < 0, 
respectively. In the limiting case a = 0, we have By 
exp[ix(1 +t?)!/2] = E_. The derivatives of the phase- 


function 6(t;x,a) with respect to t are given by 


O'(t;a) = —(@—t+2at®)/(1 + t2)!/2, (11a) 
A" (tsa) = (1—3at —2at3)/(1 + t2)3/2, (11b) 
A" (ta) = —3(at+n/(1+t2)>”2. (1c) 


The phase 6(t;a) and its first and second derivatives 


@'(t;a) and 6’ (t;a) are depicted in figures 3a, b, c, 
respectively. 


a = -1 -6 -4-2-1 a =0 06 
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Fig. 3a — The Phase Function @(t;a) for 0 < t < 6 and 
Several Values of a 
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Fig. 3b — The Function @'(t;a) for 0 < t < 6 and 
Several Values of a 


Fig. 3c — The Function 6’'(t;a) for 0 <t < 3 and 
Several Values of a 


It may be seen from figure 3b that the first 
derivative @’ does not vanish except if 0 < a < 1/23/27, 
for which @’ vanishe, for 2 values of t. Equation (ila) 
shows that the 2 points where @' =0, that is where the 
phase @ is stationary, are given by 
t,(@) = [1+(1-8a%)!/2}/4a. (12) 
We have 0 €t_ <t, € ©, witht_= 0 and t, = © for 
a = 0, andt_= 1/22 = t, fora = 1/2) The2 
points of stationary phase t_ and t, are apparent on the 
left side of figure 1 for a = .1, .2 and 1/27/. Figure 3c 
shows that the second derivative @'’ vanishes for one 


value of t, say to» if a > 0. Equation (11b) yields 


toa) = (F-1/1)/2'/7, (13) 
where [ is defined.as 
P= {(1+(14+2a2)!/472!/20}!3 (13a) 


The value ty for which @” = 0, and for which @’ reaches 
its maximum as figure 3b indicates, is a decreasing 
function of a. We have ty = © for a = 0, ty = 1/2!” 
for a = 1/23/2, and to = 0 for a = ©. Furthermore, 
figure 4, where the functions t)(a), t_(@) and t , (a) are 
depicted, shows that we have 0 <t_ <t) <t, < © for 
0<a< 1/2)” 

4 


1 | 
t+) 0.1 0.2 0.3 0.4 0.5 


a—> 


Fig. 4 — The Functions to(a), t, (a) and t_(a) 


In the vicinity of the points of stationary phase t , 
the first derivative of the phase-function may be 
approximated by the two-term Taylor series 
6'(t;a) ~ (t-t,)O% + (t-t,)°%'/2, (14) 
where the function t, (a) is defined by equation (12), and 
oe". and ie represent the values of the functions 0” (t;a) 
and 6’ (t;a) fort = t, (a). By using equation (12) in 
equations (11b) and (llc) we may obtain 
0", = +23/2a(1 —8a%)!/2/[1 + 4a? +(1—8a7)!/7]!/2, (15a) 
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OY =>-2!/296a4/[1 + 4a? + (1 — Ba%)!/7)3/. (15b) 


In the neighborhood of the point tp, where 6"= 0, we 
have 

O'(t;a) v Oy + (t—t, OG" /2, (16) 
where the function t)(a) is defined by equation (13), and 
8 and 5’ represent the values of the functions 6'(t;a) 
and @'’(t;a) defined by equations (lla) and (llc) for t = 
g(a). The functions t)(a), t , (a) and t_(a), @4(a) and 
89’ (a), O(a) and ©"(a), 8''(a) and @"" (a) are depicted 
in figure 5. The signs of the functions @/, ie oe" and 
© are readily apparent from this figure. The Taylor- 
series approximation (14) and (16) are useful for devising 
an efficient numerical method for evaluating the integrals 
(1), (a,b) and (5a,b), as will be examined in detail 


elsewhere. 


3 


Fig. 5 — The Functions to(a), t (a), Og(a), Og'’(a), 
O'/(a) and Q!/’(a) for 0 < a < 0.5 


Asymptotic approximations valid for z = 0 and x 
— —oo will now be obtained for the integrals (1), (3a,b) 
and (Sa,b). These five integrals may be expressed in the 
form 
7,(x,a) = Im(et yy +cp wy), (17) 
where 0 < k < 4, c* are constants, and w,* are the 


integrals defined as 


via) = fo expiss , (ie IKOa,(0 dt; (18) 
in this expression, tie phase 6, (t;a) is defined by 
6 (tia) = (1F at +t)!/2 with a > 0; (19) 


finally, the constants ce in equation (17) and the 
functions a,(t) in ecuation (18) are defined in the 


following table: 

ko ay yy 

0 ¢ 1 1 1 (20a) 
Oe (CREWE fi i (20b) 
Zia Ae lie t2 ye ~i (20c) 
ae 1+t? =] =f (20d) 
AO Qe ty sls ol (20e) 


In the particular case a = 0, that is on the track of 
the ship, equation (18) yields y* = yw, = yw, , with 
w,(%,0) = I explix(! + 12)!/2] K(pa,(t) dt. 

The major contribution to this integral in the limit x > 
—oo stems from the point of stationary phase at the origin 
t=0. Specifically, we may obtain 

(n/2)!/2( —x)!/24, (x,0) © 

Im c,Ky exp[i(x—1/4)] as x > —=, (21) 
with Co = 1,¢, = i, c,; =—1, (21a,b,c) 
and c, = 0 = c,; we thus have 

by —sOe— $xy for a = 0, (22a,b) 
in accordance with the symmetry of the wave pattern 
about the axis a = 0. In equation (21), Kg represents the 
value of the function K(t) at the origin t = 0, that is we 
have Ky = K(0). 

For 0 < a < 1/23”, the phase 6, (t;a) is stationary, 
that is 6’, = 0, at the two distinct points t_(a) and t, (a) 
defined by equation (12), whereas 0’ (t;a) > 0 for t > 0, 
as may be seen from figure 3b. Equations (17) and (18) 
then yield 
nd, (x,a) ¥ Im cf pt as x > —© withO<a< 1/23/2. 


The contribution of the two points of stationary phase t , 
can be evaluated by using the method of stationary phase, 
with the result 

(n/2)!/2(—x)!24, (xa) © Im o (Ay E_ +A;E,) 

as x > —0, with 0 < a < 1/23/2 and (23) 
(23a,b,c,d,e) 
furthermore, A and E, are the amplitude and 


Co = 1,¢, =1 = c, ¢; = 1 = c; 


exponential functions defined as 

M2 = Tk LEROY (24a) 
E a = exp[i(xO ¥ +n/4)), (24b) 
where K , , az, 9, and ©", are defined as 


K, = K(@), a: = at); (25a,b) 


©, = Ot,;2), ©% = 6", ;2), (25c,d) 
and t, is the function of a given by equation (12), that is 


we have 

t, = [1+(1-8a?)!/2}/4a. (26a) 
We may then obtain 

1+t,? = [1+4a?+(1—8a7)!/2}/80?, (26b) 


The expressions for the terms at may readily be obtained 
from equations (20a-e), (25b) and (26a,b). Equations (25c), 
(19) and (26a,b) yield 
O, = (3 +(1—8a%)!/2I1 + 4a? + (1 — 8a2)!/2]!/2/21/28 a. (26c) 
Finally, 0”, is given by equation (15a), that is we have 
FO" = 23/2a(1 —80%)!/2/[1 + 4a? + (1—Ba?)!/2)!/2_ (26d) 
It may be shown from equations (26a,c,d) and 
verified from figures 4, 3a, 3c and 5 that we have t_ = 0 
and @ = 1 = ©” in the limit a = O. Furthermore, 
equations (25b) and (20a-e) show that we have a, = 1 for 
k = 0, 1 and 3, and a, = 0 for k = 2 and 4 in the limit 
a = 0. The asymptotic approximation (23) for $,(x,a) 
therefore becomes identical to the asymptotic 
approximation given by equations (21), (21a,b,c) and 
(22a,b) for $,(x,0) in the limit a = 0, if the contribution 
of the second point of stationary phase t, = © is null, 
that is if A,’ = 0 for a = 0 andt, = ©. In other words, 
the asymptotic approximation (23) for $,(x,@) is uniformly 
valid in the vicinity of the track of the ship a = 0 if At = 
0 for a = 0. Equations (26a) and (26d) yield t, ~ 1/2a 
and —0", ~ 2a as a ~ 0. We then have —O”, ~ I/t, 
as a > 0, and the condition for the asymptotic 
approximation for $,.(x,a) to be uniformly valid in the 


limit a = 0 takes the form 


t!/2K(t)a,(t) > 0 as t > ~. (27) 
In the limit @ = 1/23/2, we have t_ = 1/2!/2 = 
t, and ©" = 0, as may be verified from equations 


(26a,d) and figure 5. Equation (24a) then shows that 
have |A;*| > © as a > 1/23/2; and the asymptotic 
approximation (23) is not valid in the vicinity of the 
boundary of the Kelvin wake. A complementary 
asymptotic approximation, expressed in terms of Airy 
functions, valid at and near the Kelvin cusp line is given 
in Ursell [11] for the particular case of a pressure point at 
the free surface. However, we are mostly interested in the 
sector 0 < a < 1/23/2, that is inside the Kelvin wake, in 
the present study. 

The far-field asymptotic approximation (23) shows 
that the wave pattern at any point (x,a), with x < —1 
and 0 < a < 1/23/2, consists in two elementary plane 
progressive waves. Specifically, equations (4), (20b), (23) 


and (24b) show that the free-surface elevation within the 
Kelvin wake at a sufficiently-large distance behind the ship 
is given by 
(n/2)!/2(-x)!/2e(x,a) ~~ Re[Ajexplip_) + Aj exptie , )), (28) 
where the phases ~_, of the two waves are given by 
@, =x@, +1/4. Curves along which the phases @ , or 
@_ aré constant then are defined by the equation 
xO a +n/4 = constant. This relation and equation (8) then 
yield the following parametric equations for the curves 
along which the phase is equal to —2nn: 
—x, = (2nn+7n/4)/O,, (29a) 
y, = Q@nn+n/4)a/O,, (29b) 
where 0 < a < 1/23/2. The ten constant-phase curves 
corresponding to 1 < n < 10 are depicted in figure 6. The 
“transverse”’ and “‘divergent’’ waves in this classical 
representation of the Kelvin wake correspond to the waves 
A, exp(i@ _) and AG exp(i@ , ), respectively, in equation 
(28). 


Fig. 6 — The Classical Kelvin Ship Wave Pattern 


The wavenumber corresponding to the wave with 
phase ~ , is given by V@_,. The corresponding 
wavelength, say A, and direction of propagation with 
Tespect to the track of the ship, say B.,, then are given by 
A, = 2n/|vVp,| and B, = tan~ pt /@e), where @~,* and 
Py represent the x- and y-derivatives of @, and |Vq ,| 
= ez + @; )7]!/7. The relation p, = xO, +1/4 and 
equation (26c) then yield 
A, = 2! 16na?/[3 + (1-8a?)!/7} [1-4a? + (1-8a)!/?}"”?, (30a) 
B, = sirt {[1-4o?+ (1-8a2)!/2}/[14+-4a2+ (1-8a2)!/)}!/2 (30b) 
Equations (30a,b) show that we have A_ = 2n, B_ = 0 
and A, = 0, 8B, = 1/2 in the limit a = 0, and A_ = 
4n/3 = 1, B_ = sin~1(1/3”7) = B, for a = 1/23”. 
More precisely, we have 
2n2\_ 2 4n/3 21, 20 and (31a) 
0<B_ <sin—'(1/3'7) <p, < 1/2, (31b) 
as may be seen from figure 7 where the functions 
A (a)/2n and 26 , (a)/n are depicted. 


0 0.1 0.2 0.3 1/2V2 


a— 


Fig. 7 — The Wavelengths A, (a) and Propagation Angles 
B.(a) of the Transverse and Divergent Waves in the 
Kelvin Wake 


Equation (28) shows that the amplitudes of the 
transverse and divergent waves in the Kelvin wake are 
asymptotically given by (2/n)!/2|A#* \/(—x)!?2 as x > -0. 
The steepnesses, say s,, of these waves then are given by 
8, = (2/n)!|A#|/(—x)!2, - Equations (24a), (25a,b) 
and (20b) then yield 
(-x)!7s, ~ 0, (a)|K(t,)| as x > —~, (32) 


where o , (a) is defined as (2/n)!/2(144?, )!/27( Ole) Cane 
Equations (26b,d) and (30a) then yield 
0, = B¥(1—8a2)!/I[1 — 4a + (1 —8a7)!/2)1/? 
[1 + 4a2 + (1 — 8a)!/2}3/4/ 
64n(2!/2n)!/297/2(1 — 8a2y!/4. (33) 
Equations (12) and (33) yield t_ = 0 and o_ = 
1/n(2n)!/2 in the limit «a = 0, for which we have A_ = 2n 
and B_ = 0 as was noted previously. The steepness of the 
transverse wave at a point (x,0) on the track of the ship 
then is given by 
s_(x,0) © |K(0)|/n(2n)!/2(—x)!/? as x > —20. (34) 
Equations (12) and (33) also yield t, ~ 1/2a and 
G9 1/16n3/2a’/2 in the limit a > 0, for which we have 
A,> 0 and B,> n/2. We thus have 0, » t!/2/n(2n)\/? as 
a — 0, and equation (33) shows that the steepness of the 
divergent wave at a point (x,a) in the vicinity of the track 
of the ship is given by 
5, (x,a) © t7/21K(t , )|/n(2n)!/2(—x)!/? 
asa > 0, witht, 1/2a. (35) 
The steepness s , (x,a) then becomes unbounded as 
a Oif 
t’/2|K(t)| > © as t + ©. (36) 


Equations (4), (20b) and (27) show that the asymptotic 
expansion (23) for the free-surface elevation is valid in the 


vicinity of the track of the ship if 


3/2K(t) > O ast > ©. (37) 
Let us assume that we have 

|K()| ~ 1/t# as t > ©. (38) 
Both conditions (37) and (36) are then satisfied if 

SWIPE NSU P2 (39) 


Condition (37) means that the amplitudes of the 
divergent waves in the Kelvin wake vanish as a > 0, that 
is as the track of the ship is approached, whereas 
condition (36) means that the waves become infinitely 
steep; this is’ theoretically possible because the wavelengths 
of the divergent waves vanish as a > 0, as is indicated in 
equation (3la) and figure 7. However, infinitely-steep 
water waves cannot exist in reality; indeed, there exists a 
theoretical upper bound for the steepness of water waves 
in deep water which is approximately equal to 1/7. 
Condition (39) therefore suggests that no divergent waves 
can exist within a certain domain in the vicinity of the 
track of the ship, and that the Kelvin wake contains three 
distinct regions: (i) an inner region adjacent to the track 
of the ship where only transverse waves can exist as. was 
just noted, (ii) an outer region where both transverse and 
divergent waves are present, and (ili) a region at the 
boundary between the inner and outer regions where short 
steep divergent waves, as well as transverse waves, can be 
found. It must be kept in mind, however, that these 
conclusions regarding the Kelvin wake are based on 
condition (39), which was obtained on the basis of an 
analysis in which surface-tension and nonlinearities are 
ignored. Inasmuch as this linear no-surface-tension 
analysis predicts short steep waves, both surface-tension 
and nonlinear effects are liable to be significant, and these 
effects should therefore be included in a more realistic 
analysis. In particular, it is evident from Lamb [9, 
pp. 468-470] and Wehausen and Laitone [10, pp. 636-637] 
that the system of divergent waves in the immediate 
vicinity of the track of the ship may be profoundly 


affected by surface tension. 


Yh y = £(1-»)(tan B+ z tan y) 


Fig. 8 — Waterlines and Framelines of the Simple Ship 
Bow Form Considered for Numerical Applications 


A SIMPLE TEST CASE: 
THE FAR-FIELD WAVE-AMPLITUDE FUNCTION 
The foregoing theoretical results are investigated 
numerically for the simple semi-infinite ship form studied 
previously by Scragg [5]. This ship form consists of a bow 
region, with length L, followed by a parallel body, with 
invariant framelines, extending to infinity downstream. All 
framelines, both in the bow region and downstream from 
it, are trapezoidal in shape with constant draft D. The 
waterlines are rectilinear. More precisely, the hull form is 
defined by the equations 
y = +(tanfB+z tany)(1—x) 
for0 <x <1and02z2 —d, (40a) 
y = +(tanfB+z tany) forx <0 and0 2z2 —d, (40b) 
where x,y,z and d are nondimensional in terms of the 
length of the bow region, that is we have x= X/L and 
d=D/L. Equations (40a,b) require that the condition 
tanB > d tany (41) 
be satisfied. Equation (40a) shows that the entrance angle 
at the bow (x = 1, z = 0) is equal to 28, and it may be 
seen from equations (40a,b) that y represents the flare 
angle for x < 0. The four waterlines corresponding to 
z = 0, —d/3, —2d/3, —d and the five framelines 
corresponding to x = 1, 0.75, 0.5, 0.25, 0 are depicted in 
figure 8 for d = 0.1, B = 12° and y = 45°. The 


notation 


Bo = tanB, yg = tany (42a,b) 
will be used for shortness hereafter. 

The far-field wave-amplitude function K(t) for the 
foregoing semi-infinite hull has been evaluated for two 
simple approximations defined explicitly in terms of the 
hull shape and the Froude number, namely the Michell 
thin-ship approximation and the zeroth-order slender-ship 
approximation [7], for which the function K(t) is denoted 
K\( and Kj(t), respectively. The Michell thin-ship 
approximation is given by the product of two single 


integrals, as follows: 


1 
ie = fe exp[—iv2(1 +t2)!/2x] dx 


Wh 0 2, 2 
haps exp[v*(1 +t°)z] (8) + yz) dz. 


These integrals can be evaluated analytically, with the 
result 
Ky) = 41Bp — (By — rode — F*yo(1 —e)/(1 + 1?)] 

sin{v2(1 + t?)!/2/2] exp[—iv2(1 + t2)!/2/2]/(1 + t2)3/2, (43) 


where the term e is defined as 


e = exp[—vd(1+t2)]. (44) 
Equations (43) and (44) yield 
[Ky (0| ~~ 4Bo|sin@7t/2)\/t3 as t > -. (45) 


The zeroth-order slender-ship approximation may be 
expressed in the form 
K,(t) = Kg ()+Ko @), (46) 
where Kg (t) is given by ‘the sum of a double integral over 
the hull surface and a single integral along the top 


waterline, as follows: 


4 f/f! yf 0 Din se1D 
Ky = vf, dx ar dz exp[v*(1 + t*)z] 

exp[—iv?(1 + t7)!/2(x + yt)](Bp + ¥o2) 

fi! dx exp[—iv2(1 + t?)!/2(x + yt)] 

Ba/ [1 + Bh + ¥601 — X71. (47) 
The integration in these integrals is carried over the 
positive half of the hull surface. Equations (40a) and 
(42a,b) then yield y = (By +yz)(1 —x). We then have 
xtyt = 1—(1 +fot)(1 —x) + yot(1 — x)z. 
By using this relation into equation (47) we may then 


obtain 
expfiv(1 +t?)!/7]K F(t) = vf} exp[iv?(1 + t?)!/2(1 +fot)x] 


“ [Bolg —yolj* — B3/(1 + B2 + vex?)] dx, (48) 


where I with n = 0 and |, are the integrals defined as 


(-)t = v7? fs a exp[v>(1 +t?(1 + iygux)z}z" dz, 
with u defined as 
u = t/(1+t?)!?, (49) 
The integrals Ie may be evaluated analytically with 
the result 
Ig = (1-ee,)/(1+t?)(1 Fiygux), 
I* = [F°(—ee,)/(1 +t?)(1 ¥ iygux)— dee , ] 
/(1 +t)(1 Fiygux), 
where e is the exponential function given by equation (44) 
and e, is the exponential function defined as 
e, = exp[+iv7dygt(1 +t?)!/2x). (50) 
We then have 
Bolg —Yolj —Bg/(1 + BG + vox") = 
A, —ee ,B, /(1 +t?)(1 ¥ iygux), (51) 
where the terms A, and B, are given by 
A, = Bo/(1 +t)(1F iygux) — B3/(1 + B3 + ax") 
—yF*/(1 +1?)7(1F iygux)?, 
B, = By—¥od—yoF?/(1 + t?)(1 + iygux). 
It may be verified that we have 
A, = Bol(1 —Bet?)/(1 +t?) + x7B2¥2/(1 + BS + ya) + Bd) 
yg Coes (52a) 
B, /(1+t?)(1F iygux) = Bo/(1+t*)—yoD,, (52b) 
where the terms C, and D, are given by 
C= [Baya/(1 + Bay + BS + volx7(d — x°)/( + BS + 6x7) 
— [Boy ot2x? + F2(1 + t? — yot?x)/(1 + t? + y@t?x?) 
+ it(1 + t?)!/2{p5—2F?yq/(1 +t? + yat?x?)}x] 
JA +t?\(1 +t? + y2t?x?), (53a) 
(1+?) +t? +yAt7x7)D, = Boyot?x? + d(1 +1) 
+P(1 +0? —ydt?xy/(1 +t? + yet?x?) 
F it(1+t?)!/2[B)—yod — 2F*y9/(1 + t?+ yot?x?)]x. (536) 
Equation (50) yields 
exp[iv2(1 + 17)! + Bot)xle, = 
expfiv(1 + t?)!/(1 = ByDxl, (55) 


where f3, is defined as 


By = Bo—dy¥o = tanf—dtany, (56) 
as may be obtained from equations (42a,b). By using 
equations (55), (51), (52a,b) and (44) into equation (48) we 
may obtain 

expfiv2(1 +17)! 7K (0) = 

Bol(1 — Bt?)C1 +t?) ~ "1G (Bo) + Bayo! + Bo + YO) '1;* (sBo) 

— (1+ Bdyexp{ —v2d(1 + t?)}(1 + t?)~ Ne (5B IZ + 82) 
+yolJ , (tsBy.C .) + exp{ —v7d(1+ t?)}J, (6By,D,)1, (57) 
where the functions I (t;B) and J, (t;B,A), or more 


precisely I*(t;8,v) and J, (t;8,A,v), are defined in terms of 
the integrals 


1 (t;8) = ye exp[iv2(1 +t2)!/2(1 + pt)x]}x2" dx, (58a) 


J,.(tsB,A) = v? yf exp[iv2(1 +t2)!/2(1 = pt)xJA(x) dx. (58b) 


The integrals I* (t;8) can be evaluated analytically, 
with the result 
Ig (t;B) i(I-E,)/o,, (59a) 
IF (t;B) —if(i+ 2iF*/o,, JES + 2F4(1-E - \/o?, \/o., , (59b) 
where o, and E, are defined as 
o, = (1+t?)!(1 pt), (59c,d) 
Expressions (59a,b) for Ig (t;B) and IG (t;B) are not valid 
1/B, for which 


E += exp(iv2o P): 


in the special case when we have t = 
expressions (59a,b) become 
If (1/B;8) = v?, 1; 1/B;B) = v2/3. (59e,f) 
The integrals J, defined by equation (58b), where 
the amplitude function A(x) takes the form of C.. or D, 
specified in equations (53a,b), cannot be evaluated 
analytically. These integrals were then evaluated 
numerically by dividing the integration range 0 < x < 1 
into N segments of equal length and using piecewise 
quadratic approximations for the amplitude function A(x) 


within each segment. In this manner, we may obtain 


N 
oJ, (t;B,A) Ze) MWA ej.) 


+ (NF*/o , (e, — 1A; , 1 — Aj) + QNF*/o ,.) 
(ler. +i(2NF*/o , le, - D}(A\+A;, 1-24; 4 1/2)) (60) 
where o, is given by equation (59c), , is defined as 


(60a) 


+ 


exp(ivo , /N), 
A 


5a = 


and vp A 


AG i represent the values of the 


G-D/N, % 44 


= (j—1/2)/N, respectively. Expression 


j+1/2 
amplitude function A(x) at the points ie = 
= j/N, and Xi 41/2 
(60) for J, (tB,A) is not valid in the special case when t 


= 1/8, for which we have 


N 
6NJ , (1/B;B,A) © 2 A,+A,,,+4A (60b) 
j= 


In summary, the zeroth-order slender-ship 


j+1/2° 


approximation K,(t) is determined by equations (46), (57), 
(59a-f), (60), (60a,b), (42a,b) and (56). In the limiting case 
Yo = 0, these equations yield 

Ko(t) = 4Bo(1 +12) 3/7((1 — pat?) + pd) — 


—exp{ —v7d(1 + t7)}](R —i1)/(1 — p2t?), (61) 
where the terms R and I are given by 
R = sin{v2(1 + t2)!/2/2] cos[v2(1 + t2)!/2/2] 
— foot sin[v?Bot(1 + t?)!/2/2] cosfv°Bot(1 +t7)!/2/2], (61a) 
I = sin?[v?(1 + t?)!/2/2] —sin?[v7p,t(1 + t7)!/2/2]. (61b) 
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In the thin-ship limit By ~ 0, equations (61) and (61a,b) 
yield 

Ko(t) © 48o(1 +t?) 3/2[1 — exp{ — v7d(1 + t?)}). 

sinfv(1 + t?)!/2/2) exp[—iv2(1 +1?)!/2/2] as By > 0. 
The thin-ship limit (62) of the zeroth-order slender-ship 


(62) 


approximation K,(t) may be seen to be identical to the 
Michell thin-ship approximation (43) in the particular case 
Yo = 0. However, the thin-ship limit (62) of the slender- 
ship approximation (61) is not uniformly valid in the limit 
t — ©; indeed, equations (61) and (6la,b) yield equation 
(62) in the limit Bot > 0. Equation (62) yields 


|Ko(t)| ~ 4Bolsin(v2t/2)|/t? for 1K t K 1/fp. (63) 
More generally, equation (61) yields 
[Ko(t)| ~ 4B9|R —il|/(1 + B6)t3 as t + ©. (64) 


Equations (61a,b) show that |R—il| ~ Bot|sin(v2B,t?)|/2 if 

Bot >> 1. We then have 

|Ko(t)| ~~ 2B2|sin(v2Bot?)|/(1 + Bayt? for t > 1/By. (65) 
In the limit t~°°, equations (57) and (59a-d) yield 

expliv*t)Kg (t) © Bol(l — BGt7VIg (tsBo) 

+ yQ(1 + BS +76) 'BGt7T;! (tsBg)I/(1 + Bot as t > ©, 

where we have 

Ig (tsBo) © if — exp{ivt( + Bot)}]/(1 ¥ Bott as t > %, 

I}* (tsBo) ‘Y —i exp{iv7t(I F Bot)}/(1 + Bott as t > &. 

We then have 


|Ko(t)| © 2Bo|N|/(1 + B9)|1 — Bat |t? as t + -, (66) 
where the term N is given by 
N = (1—ft?)sin(v2t) +i cos(v’t)) 

— (1-07 B?t?)[Bot sin(vBot?) +i cos(vBot”)], (66a) 
with o? defined as 
o? = (1+d)/(1 +2 + yA). (66b) 
We may then obtain 
[K(t)| ~~ 4Bo|sin(v2t/2)|/(1 + BA) for 
1<t <F/G,)!, (67a) 


[Ko(t)| ~ 2A2|sin(v2Bot?)|/(1 + BS + yo)t? for t >> 1/Bp. (67b) 
Equations (67a,b) are identical to equations (63) and (65) 
in the limits By > 0 and yg > 0, respectively. More 
generally, equation (66a) yields the following upper bound 
for the term |N| in equation (66) 
|N|2 < (1 —pat?)? + (1+ B2t2)(1 — 08217)? 
+ 2(1 + Bgt)?(1 + oBot)|(1 — Bot)(1 — oBgt)|. 
Equations (67a,b) thus show that in the limit t > © 
we have |Ko(t)| ~ 1/3 forl<«Kt<« F/(Bp)'”, that is for 
moderately large values of t, and |Ko(t)| ~ 1/? for 1/Bo 


(68) 


< t, that is for very large values of t. The asymptotic 
approximations (67a,b) and (45) show that the Michell 
approximation K,,(t) corresponds to the thin-ship limit Bo 
< 1 of the zeroth-order slender-ship approximation K,(t), 
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Fig. 9a — Real and Imaginary Parts of the Functions Fig. 9b — Real and Imaginary Parts of the Functions 
Ky() and Ko(t) for a Simple Ship Bow Form with Kyy(t) and Ko(t) for a Simple Ship Bow Form with 
B = 12° (Top) and 3° (Bottom), y = 0 and F = 0.5 B = 1° (Top) and 20’ (Bottom), y = 0 and F = 0.5 
and that this thin-ship limit is not uniformly valid in the for values of B equal to 1° and 20’, which are quite 
limit t > ©. More generally, the limiting processes By > 0 small, and for values of t that are much smaller than 
and t > © cannot be interchanged. 1/Bo. 
Figures 9a,b depict the real and/or imaginary parts The top and bottom parts of figure 10 depict the 
of the functions (1 +17)!7K (0) and (1 +10?)1/2K (t) for the real and imaginary parts of the functions (1 +07)1K (t) 
simple bow shape defined by equations (40a,b), where the and (1 +t2)!/ *Ko(t), respectively, for the previously- 
nondimensional draft d and the maximum flare angle y considered simple bow shape with d = 0.1, B = 12°, F 
are taken equal to d = 0.1 and y = O, and four values = 0.5 and for two values of the flare angle, namely for y 
of the half-entrance-angle f are considered, namely B = = 0 and 45°. The top part of the figure shows that 
12°, 3°, 1° and 20’. The Froude number based on the differences between the curves corresponding to y = 0 
length of the bow region is taken equal to F = 0.5 in the and y = 45° are faily small, and are appreciable only for 
numerical results presented in figures 9a,b and in figures small values of t, for the Michell thin-ship approximation 
10 and 11 considered further on. The values of 1/f) = Ky. In particular, the asymptotic approximation 
1/tanfB corresponding to the values of f equal to 12°, 3°, Ky(t) © 4B sin(v2t/2) exp(—iv2t/2)/t3 ast > ©, 
1° and 20’ are approximately equal to 4.7, 19, 57 and which may be obtained from equations (43) and (44), is 
172, respectively. The functions (1 +t7)!/2K(t) and independent of y. The bottom part of figure 10 shows 
(a +1?)!/2K (t) are depicted for 0 < t < 10 in figures 9a,b. differences between the curves corresponding to y = 0 
Differences between the approximations K,, and Ko, and y = 45° for the slender-ship approximation Ky that 
especially their imaginary parts, can be seen to be are significantly larger than those for the Michell 
substantial in figure 9a corresponding to B=12° and 3°. approximation K,,, especially for intermediate values of t 
Figure 9b shows that differences between the imaginary in the vicinity of t = 1/tanB ~ 4.7. The flare angle y 
parts of the functions K,, and Kg remain appreciable even thus has a pronounced effect upon the behavior of the 


iit 


f B = 12°, F = 0.6 
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Fig. 10 — Real and Imaginary Parts of the Functions 
Kyy(t) and Ko(t) for a Simple Ship Bow Form with 
y = O and 45°, B = 12° and F = 0.5 


slender-ship approximation K,(t) for values of t in the 
neighborhood of 1/tanf. The influence of y upon |Ko(t)| 
for large values of t is explicitly indicated by the 
asymptotic approximations (67a,b). 

Finally, figure 11 depicts the real and imaginary 
parts of the functions (1 +1?7)'/K(t) and (1 +1?)!/2K (t) 
for the simple bow shape depicted in figure 8, for which 
we have d = 0.1, B = 12° and y = 45°, at a value of 
the Froude number F equal to 0.5. Differences between 
the approximations K,, and Ky may be seen to be quite 
substantial. In particular, the function (1 +12)!/ *Ko(t) has 


a peak at t = 1/tanf = 1/fp. 
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B = 12°, y = 45°, F=0.5 


Fig. 11 — Real (Top) and Imaginary (Bottom) Parts of 
the Functions K,4(t) and Ko(t) for a Simple Ship Bow 
Form with B = 12°, y = 45° and F = 0.5 


A SIMPLE TEST CASE: 
THE FAR-FIELD KELVIN WAKE 


The expressions for the far-field wave-amplitude 
function K(t) obtained in the foregoing section for a 
simple bow shape and for two simple approximations to 
the function K(t), namely the Michell thin-ship 
approximation K,, and the zeroth-order slender-ship 
approximation Kp, may now be used into the previously 
determined asymptotic approximations for the far-field 
Kelvin wake. 

Far behind the ship, that is for x > —°, the 
functions $4 (x), where 0 < k < 4, defined by equations 
(17), (18), (19) and (20a-e) are given by the asymptotic 
approximation (23). The real and imaginary parts of the 
amplitude functions A, (@) and Ay (a) in this asymptotic 
approximation, given by equation (24a), are depicted in 
figures 12a and b, 13a and b and 14 for 0 < a < 1/23/. 
More precisely, figures 12a,b and 13a,b represent the 
amplitude functions A, and A,’ for k = 0,1,2,3 
associated with the potential ¢ and its derivatives $,, ty 
and 9x49 and correspond to the approximations K,, and 


Kp, respectively. Figure 14 depicts the amplitude functions 
A, and Aj corresponding to txy for the approximations 
Ku and Kp. 
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Fig. 12a — Real (Top) and Imaginary (Bottom) Parts of 
the Functions K af (+04)! for k = 0 (Left) and 1 
(Right), K = Ky (Michell Thin-Ship Approximation) and 
a Simple Ship Bow Form with B = 12°, y = 45° and 
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Fig. 12b — Real (Top) and Imaginary (Bottom) Parts of 
the Functions K , ag /(+04,)!/2 for k = 2 (Left) and 3 
(Right), K = Ky, (Michell Thin-Ship Approximation) and 
a Simple Ship Bow Form with B = 12°, y = 45° and 
|? = (0.5 

It may be seen from figures 12a and b, 13a and b 
and 14 that the amplitude functions Ay (a) and AX (a) 
become unbounded in the limit a > 1/27/7. This 
singularity at a = 1/23/2 stems from the fact that the 
asymptotic approximation (23) is not uniformly valid in 
the limit a ~ 1/23/2, as was already noted. The amplitude 
functions Ay (a) corresponding to the system of divergent 
waves in the Kelvin wake also become unbounded in the 


limit a > 0 in figures 13b and 14 corresponding to k = 2 
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Fig. 13a — Real (Top) and Imaginary (Bottom) Parts of 
the Functions K , ag /( +0 {)!/2 for k = 0 (Left) and 1 
(Right), K = Ko (Zeroth-Order Slender-Ship 
Approximation) and a Simple Ship Bow Form with 
B = 12°, y = 45° and F = 0.5 
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Fig. 13b — Real (Top) and Imaginary (Bottom) Parts of 
the Functions K , a /(+ 9;)!/2 for k = 2 (Left) and 3 
(Right), K = Ko (Zeroth-Order Slender-Ship 
Approximation) and a Simple Ship Bow Form with 
Bp = 12°, y = 45° and F = 0.5 
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and 3 for the slender-ship approximation Ky and to k = 4 
for both approximations K,, and Ko, respectively. This 
singularity at a = 0 illustrates the previously-noted 
conclusion that the asymptotic approximation (23) is not 
uniformly valid in the vicinity of the track of the ship a 
= 0 if condition (27) is not satisfied. Equations (20a-e) 
show that we have a) © 1, a, ~ t, a) © (7, a, © t? and 
a, © t? as t > ©, and equations (45) and (67b) yield 
|Kyl ~ 1/3 and |Kyl ~ 1/1? as t > ©, respectively. 
Condition (27) therefore is not satisfied for k 2 4 and k 
> 2 for the approximations K,, and Ko, respectively. 
Condition (27) however is satisfied for k = 0 and 1, 


corresponding to the potential $ and the free-surface 
elevation $,, for both the approximations Ky, and Kp. 


Figures 12a and 13a show that the amplitude 
functions Ay and A; , corresponding to the transverse 
waves in the Kelvin wake, generally are larger in 
magnitude than the amplitude functions Ag and At 
corresponding to the divergent waves, whereas the reverse 
may generally be seen to hold in figures 12b and 13b for 
the amplitude functions A;* and Af. The relative 
importance of the divergent waves with respect to the 
transverse waves thus increases with k. Indeed, the 
transverse-wave amplitude function A, (a) is hardly visible 
on the scale of the divergent-wave-amplitude function 
Aj (a) used in figure 14. 
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Fig. 14 — Real (Top) and Imaginary (Bottom) Parts of 
the Functions K , a¢ /( #01/)!/2 for Kyy (Left) and Kg 
(Right) and a Simple Ship Bow Form with B = 12°, 
y = 45° and F = 0.5 
The divergent-wave-amplitude functions Ay (a) 

associated with the approximation Ky are most notably 
different from the corresponding functions Aft (@) 
associated with the Michell approximation K,, for values 
of a in the vicinity of a = 0, as was already noted, and 
tanB/(2 + tan’). In the vicinity of this value of a, 
the divergent-wave-amplitude functions Ay (a) associated 


of a = 


with the approximation K, exhibit a peak, which is quite 
pronounced for k 2 1. The foregoing particular value of 
a corresponds to the special case when the point of 
stationary phase t, , defined by equation (12), is equal to 
the value 1/8) = 1/tanf for which the function Ko(t) 
displays a peak, as may be seen from figures 10 and 11. 
Figure 15 depicts the amplitude functions 
(2/n)'/2AP and (2/n)!/2A* that are associated with the 


free-surface elevation far behind the ship, as is specifically 
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Fig. 15 — Amplitude of the Transverse (—) and 
Divergent (+) Waves in the Kelvin Wake for a Simple 
Ship Bow Form, with B = 12°, y = 45° and F = 0.5, in 
the Michell Thin-Ship Approximation (Left) and the 
Zeroth-Order Slender-Ship Approximation (Right) 
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indicated in equation (28), for the previously-considered 
ship bow shape. It may be seen that the amplitude a , (a) 
of the divergent waves in the Kelvin wake vanishes as a > 
0, that is at the track of the ship, and is generally smaller 
than the amplitude a_ (qa) of the transverse waves; this is 
especially true in the vicinity of the track of the ship. 
Differences between the wave-amplitude functions a , (a) 
corresponding to the approximations K,,(t) and Ko(t) are 
particularly striking for the amplitude a, (a) of the 
tanB/(2 +tan2p), 
where the function a , (a) associated with the 


divergent waves in the vicinity of a = 


approximation K)(t) exhibits a sharp peak. 

Figure 16 depicts the wave-steepness functions s_ (a) 
and s , (a), which correspond to the ratios of the wave- 
amplitude functions a_(a) and a, (a) depicted in figure 15 
over the wavelength functions A_(a) and 1, (a) defined by 
equation (30a); the steepness functions s , (a) are 
specifically defined by equations (32) and (33). The 
divergent waves in the Kelvin wake may be seen to be 
generally steeper than the transverse waves, even though 
figure 15 shows the transverse waves to be larger in 
amplitude than the divergent waves. This is especially true 
near the track of the ship where the steepness of the 
divergent waves becomes infinitely large, even though 
figure 15 shows that their amplitude vanishes as a > 0. 
The divergent waves become infinitely steep at the track 
of the ship because the wavelength A, (a) > 0 as a > 0 
and condition (36) is satisfied, for both the 
approximations K,,(t) and K,(t) as may be seen from 
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Fig. 16 — Steepness of the Transverse (—) and Divergent 
(+) Waves in the Kelvin Wake for a Simple Ship Bow 
Form, with B = 12°, y = 45° and F = 0.5, in the 
Michell Thin-Ship Approximation (Left) and the Zeroth- 
Order Slender-Ship Approximation (Right) 
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Fig. 17 — Upper Bound for the Steepness of the 
Divergent Waves in the Kelvin Wake for a Simple Ship 
Bow Form, with 8 = 12° and y = 45°, in the Michell 
Thin-Ship Approximation (Left) and the Zeroth-Order 

Slender-Ship Approximation (Right) 
equations (45) and (67b). Figure 16 also shows a sharp 
peak in the steepness of the divergent waves at the value 
of a equal to tanf/(2 + tan). 

The steepness s , (a) of the divergent waves is given 
by equations (32) and (33). An upper bound for the 
functions s , (a) may be obtained by using an upper 
bound for the function |K(t ,)| in equation (32). Equation 
(45) yields the following upper bound for the function 
|Ky(0| corresponding to the Michell thin-ship 
approximation: |Ky(o| < 4Bo/t? as t ~ 9. An upper 
bound for the function |Kp(t)| is given by equation (66), 
where the upper bound defined by equation (68) is used 
for the term |N|. These upper bounds for the functions 


115) 


|Ky4())| and |K,(t)| can be expressed in terms of a by 
using equation (12). The corresponding upper bounds for 
the steepness function s , (a) are depicted in figure 17. 
Comparison of figures 16 and 17 shows that the upper 
bound for the steepness of the divergent waves in the 
Kelvin wake depicted in figure 17 is satisfactory for all 
values of @ for the Michell approximation Ky, whereas 
that corresponding to the slender-ship approximation Ky is 
satisfactory for values of a smaller than approximately 
half the value tanB/(2+tan2f). In both cases, the upper 
bounds for the function s , (a) depicted in figure 17 are 
satisfactory for the range of small values of a for which 
the steepness s , (a) is large. It is noteworthy that these 
upper bounds for the steepness of the divergent waves are 
valid for all Froude numbers, since equations (12), (32) 
and (33), and the upper bounds for |K,,(t)| and |K)(t)| do 
not involve the Froude number. 

Figure 18 depicts the steepness function s , (a) of 
the divergent waves corresponding to the slender-ship 
approximation K,(t) for the simple bow shape considered 
previously in the two cases when the maximum flare angle 
y is taken equal to 0 and 45° (at the top and bottom 
halves of the figure, respectively) and for two values of the 
Froude number, namely for F = 0.8 and 0.3 (on the right 
and left halves of the figure, respectively). Comparison of 
the right and left halves of figure 18 shows no appreciable 
difference between the values of the steepness function 
s , (a) for the range of small values of a for which the 
steepness is large, in agreement with the previously-noted 
result that the upper bound for s , (a) depicted on the 
right side of figure 17 is independent of the Froude 
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Fig. 18 — Steepness of the Divergent Waves in the Kelvin 
Wake for a Simple Ship Bow Form with B = 12°, y = 0 
(Top) and 45° (Bottom), and F = 0.3 (Left) and 0.8 
(Right) in the Zeroth-Order Slender-Ship Approximation 
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Fig. 19 — Amplitude of the Peak Value of the Steepness 
of the Divergent Waves in the Kelvin Wake for a Simple 
Ship Bow Form with 8: = 12° and y = 45° in the 
Zeroth-Order Slender-Ship Approximation 

number. Comparison of the top and bottom halves of 
figure 18, corresponding to y = O and 45° as was already 
noted, shows appreciable differences, especially for F = 
0.3 for which the function S, (a) exhibits a very 
pronounced peak in the case y = 45°. No such peak is 
apparent in the top half of the figure for y = 0, that is 
in the case when the hull intersects the free surface 
orthogonally. 

The right half of figure 16 and the bottom half of 
figure 18 show that the magnitude of the peak in the 
steepness function s ,(a) fora = tangp/(2 + tan?) strongly 
depends on the value of the Froude number. Specifically, 
the peak is very pronounced in figure 18 for F = 0.3, 
fairly pronounced in figure 16 for F = 0.5, and almost 
nonapparent in figure 18 for F = 0.8. The magnitude of 
the peak in the steepness function, that is the value of the 
function s , (a) for a = tanf/(2+ tan’), is represented in 
figure 19 as a function of the Froude number, which is 
based on the length of the bow region. This figure shows 
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that the magnitude of the peak increases very rapidly as 
the Froude number decreases below a certain threshold — 
value in the vicinity of F = 0.3. 

Figure 20 depicts the boundary of the Kelvin wake, 
—y/x = 1/23/2 (that is, an 
angle equal to approximately 19°28’), the line a = 
tanf/(2+tan?B) ~ f/2 (that is, an angle equal to 
approximately 6°) along which the steepness of the 


which corresponds to a = 


divergent waves has a peak, and the lines along which the 
steepness of the divergent waves is equal to 1/20, 1/15, 
and 1/7 (shown as a chain line close to the track of the 
ship). The latter three lines were determined by using the 
upper bound for the steepness function s , (a) that was 
determined previously and depicted on the right half of 
figure 17. The four lines inside the Kelvin wake shown in 
figure 20 correspond to the zeroth-order slender-ship 
approximation K,(t) for the simple bow shape considered 
previously, with B = 12°, y = 45° and d = 0.1. The 
three lines along which the divergent waves are steep lie 
much closer to the track of the ship than the lines a ~ 6° 
along which the steepness of the divergent waves exhibits 
a peak. 
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Fig. 20 — The Kelvin Cusp Line, the Line Along Which 

the Amplitude of the Divergent Waves Exhibits a Peak, 

and the Three Lines Along Which the Steepness of the 

Divergent Waves is Equal to 1/20, 1/15 and 1/7 (Chain 

Line Close to the Track of the Ship) for a Simple Ship 

Bow Form, with B = 12° and y = 45°, in the Zeroth- 
Order Slender-Ship Approximation 


The three lines along which the steepness of the 
divergent waves is equal to 1/20, 1/15 and 1/7, which are 
depicted in figure 20 for 0 > x > —100 and 0 < y < 40, 
are represented again in figure 21 at a distorted scale 
where 0 > x > —300 and 0 < y < 7. The corresponding 
constant-steepness lines predicted by the Michell thin-ship 
approximation are also shown in figure 21 for 
comparison. The latter lines were determined from 
equations (32) and (33) and the upper bound |K,,(t)| < 
4Bp/t3 given by equation (45). Figure 21 shows significant 
differences between. the constant-steepness lines predicted 


Fig. 21 — Lines Along Which the Steepness of the 
Divergent Waves is Equal to 1/20, 1/15 and 1/7 for a 
Simple Ship Bow Form, with 8 = 12° and y = 45°, in 
the Michell Thin-Ship Approximation (— — —) and the 

Zeroth-Order Slender-Ship Approximation ( ) 


Fig. 22 — Lines Along Which the Steepness of the 
Divergent Waves is Equal to 1/15 for a Simple Ship Bow 
Form with B = 8°, 12°, 16° and y = 45° (Top) and 
y = 0, 25°, 45° and 8B = 12° (Bottom) in the Zeroth- 
Order Slender-Ship Approximation 


by the Michell thin-ship approximation K,,(t) and the 
zeroth-order slender-ship approximation K)(t). This figure 
strongly suggests the need for performing additional 
calculations based on a more realistic mathematical model 
than the simple thin-ship and slender-ship approximations 
used in this study. 

Finally, the effect of the entrance angle fB and of 
the flare angle y on the steepness of the divergent waves is 
illustrated in figure 22. Specifically, this figure depicts the 
lines along which the steepness of the divergent waves, as 
predicted by the slender-ship approximation Ko(t), is equal 
to 1/15 for B = 8°, 12°, 16° and y = 45°, in the top 
half of the figure, and for y = 0, 25°, 45° and B = 12° 
in the bottom half of the figure. This figure shows that 
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the short divergent waves in the Kelvin wake become 
steeper as the entrance angle f increases and as the flare 
angle y decreases. More generally, figure 22 shows that the 
short divergent waves in the Kelvin wake are strongly 
influenced by the hull shape, and it therefore suggests the 
need for performing additional calculations in which 
systematic variations in hull shape are considered. 


SUMMARY OF RESULTS AND CONCLUSIONS 

Asymptotic expressions for determining the velocity 
potential and its derivatives at a sufficiently-large distance 
behind a ship advancing at constant speed in calm water 
are given by equations (23), (23a-e), (24a,b), (25a-d), (20a- 
e) and (26a-d). The far-field asymptotic approximation 
(23) is uniformly valid in the vicinity of the track of the 
ship a = 0 if condition (27) is satisfied. For the simple 
bow shape considered in this study, condition (27) is 
satisfied for k = 0 and 1, corresponding to the potential 
$ and the free-surface elevation $,, when the far-field 
wave-amplitude function K(t) is approximated by the 
Michell thin-ship approximation K,,(t) or the zeroth-order 
slender-ship approximation K,(t). However, condition (27) 
is not satisfied for k > 4 and k > 2 for the 
approximations K,, and Kp, respectively. 

The asymptotic approximations used in this study 
provide simple explicit analytical expressions for 
determining the velocity potential and its derivatives for 
large values of — Xg/U2, that is in the far field, in terms 
of the far-field wave-amplitude function. However, for 
small and intermediate values of Xg/U2, these asymptotic 
approximations are not useful, and the integrals (18) must 
be evaluated numerically. For intermediate values of 
Xg/U2, the exponential function E, (x) +E_(t3x) in the 
integrands of the integrals (1), (3a,b) and (Sa,b) oscillates 
fairly rapidly, as may be seen from figure 1. Accurate and 
efficient integration rules suited to oscillatory integrands of 
the type depicted in figure 1 must be used. For small 
values of Xg/U2, on the other hand, the oscillations of 
the exponential function E , (t;x)+E_(t;x) are not 
significantly more rapid than the oscillations of the far- 
field wave-amplitude function K(t) which also appears in 
the integrands of the wave integrals (18), so that a 
different integration rule is required. 

The amplitude, a, , of the divergent waves in the 
Kelvin wake vanishes at the track of the ship if condition 


(37) is satisfied. However, it is well known that the 


wavelength, A, , of the divergent waves also vanishes at 
the track of the ship, as may be seen from figures 6 and 
7. Therefore, the divergent waves can theoretically become 
infinitely steep at the track of the ship. More precisely, 
the steepness, s, = a, /A,, of the divergent waves is 
unbounded at the track of the ship if condition (36) is 
satisfied. Both conditions (37) and (36) can be satisfied 
simultaneously if condition (39) is satisfied, where the far- 
field wave-amplitude function is of order 1/t4 as t > ©. 
Conditions (37) and (36), and consequently also condition 
(39), are satisfied in the cases of the thin-ship and the 
slender-ship approximations K,, and Kp for the simple 
ship form considered in the study, as may be verified 
from figures 15 and 16 where a, (a) > 0 and s, (a) > © 


asa 0. 


Infinitely-steep water waves cannot exist in reality. 
Indeed, there exists a theoretical upper bound for the 
steepness of water waves in deep water, which is 
approximately equal to 1/7. Condition (39) thus suggests 
that no divergent waves can exist within a certain region 
in the vicinity of the track of the ship, and that the Kelvin 
wake contains three distinct regions: (i) an inner region 
adjacent to the track of the ship where only transverse 
waves can exist, (ii) an outer region where both transverse 
and divergent waves are present, and (iii) an intermediate 
region at the boundary between the inner and outer 
regions where steep short divergent waves, as well as 
transverse waves, can be found. 

Surface-tension and nonlinear effects have been 
ignored in the analysis presented in this study. This linear 
no-surface-tension analysis predicts extremely short and 
steep waves in the vicinity of the track of the ship. Both 
surface-tension and nonlinear effects therefore are liable to 
be significant, and these effects should be taken into 
account. A linear analysis including surface-tension effects 
should be performed first, since it is evident from the 
results obtained in the present study and from the brief 
description of the effects of surface tension upon the 
Kelvin wake given in Lamb [9, pp. 468-470] and 
Wehausen and Laitone [10, pp. 636-637] that the system 
of divergent waves in the vicinity of the track of the ship 
is likely to be profoundly affected by surface tension. 

It was previously found by Scragg [5] that, for a 
ship bow form with a large flare angle, the zeroth-order 
slender-ship approximation K,(t) predicts a sharp peak in 
the value of the amplitude of the divergent waves at a 
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value of a equal to approximately half the entrance angle 
B. This finding of Scragg has been verified in this study, 
as may be seen from figures 15, 16 and 18. Furthermore, 
the magnitude of the steepness of the divergent waves has 
been found to increase very rapidly as the Froude number 
decreases below a certain threshold value, as is shown in 
figure 19. 

The line along which the steepness s , (a) of the 
divergent waves has a peak and the lines along which 
S, (a) takes the large values 1/7, 1/15 and 1/20 have been 
determined, for a simple ship form, on the basis of both 
the zeroth-order slender-ship approximation K,(t) and the 
Michell thin-ship approximation K,,(t). Figure 20, 
corresponding to the slender-ship approximation Ko(t), 
shows that these lines are well inside the Kelvin angle, and 
that the large-steepness lines are much closer to the track 
of the ship than the line corresponding to the peak in the 
steepness of the divergent waves. 

The lines, depicted in figures 20 and 21, along 
which the steepness of the divergent waves takes large 
constant values are independent of the value of the 
Froude number, but they strongly depend on the hull 
shape, as may be seen from figure 22 where ‘‘constant- 
steepness lines’’ corresponding to several values of the 
entrance angle f and of the flare angle y are depicted. 
This figure shows that the short divergent waves in the 
Kelvin wake become steeper as the entrance angle 
increases and/or as the flare angle decreases. 

It was found that the lines, along which the 
steepness of the divergent waves takes large constant 
values, predicted by the slender-ship approximation K,(t) 
and the thin-ship approximation K,,(t) are quite different 
from one another, as may be seen from figure 21. This 
figure therefore indicates the need for performing 
additional calculations based on a more realistic 
mathematical model than the simple thin-ship and slender- 
ship approximations used in this study. These two 
approximations correspond to simple special cases of the 
Neumann-Kelvin theory, which should then be used. In 
particular, it would be useful to determine whether this 
more realistic theory predicts that the steepness of the 
divergent waves in the Kelvin wake exhibits a peak (or 
several peaks), as was found by using the slender-ship 
approximation K,(t) for a ship bow form with large flare 
angle. Figure 21 specifically demonstrates the importance 
of obtaining accurate predictions of the far-field wave- 


amplitude function K(t) for large values of t. Indeed, the 
short divergent waves in the vicinity of the track of the 
ship, that is for small values of a, are associated with tc 
value of the function K(t) for large values of t, as may be 
seen from equations (23), (24a), (25a) and (12). Precise 
knowledge of the asymptotic behavior of the function K(t) 


as t > © therefore is critical. 
The behavior of the function K(t) as t > © has 


been determined analytically in this study for the simple 
case of the thin-ship and slender-ship approximations 
Ky and K,(0) for an idealized ship bow form. More 
precisely, the asymptotic behavior of the functions K,,(t) 
and K,(t) for the simple ship bow shape considered in this 
study is specified by equations (45) and (67a,b). These 
asymptotic approximations show that we have [Kyl = 
0(1/t3) and |Ko(t)| = 0(1/t2) as t > ©. They also show 
that the Michell thin-ship approximation K,,(t), which 
corresponds to the thin-ship limit of the slender-ship 
approximation K,(t), is not uniformly valid in the limit t 
— ©. It is possible to analytically determine the behavior 
of the far-field wave-amplitude function K(t) associated 
with the Neumann-Kelvin theory for an arbitrary ship 
form, as is shown in [12]. Such an asymptotic 
approximation for the function K(t) for large values of t 
is useful because it provides an explicit analytical 
relationship between the hull form and the Froude 
number, on one hand, and the amplitude of the short 
divergent waves in the vicinity of the track of the ship, on 


the other hand. 
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